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ABSTRACT 

Aims. Dwarf irregular and blue compact galaxies are very interesting objects since they are relatively simple and 
unevolved. We aim at deriving the formation and chemical evolution history of late-type dwarf galaxies, and compare 
it with DLA systems. 

Methods. We present new models for the chemical evolution of these galaxies by assuming different regimes of star 
formation (bursting and continuous) and different kinds of galactic winds (normal and metal-enhanced). The dark-to- 
baryonic mass ratio is assumed to be 10 in these models. The chemical evolution model follows in detail the evolution 
of He, C, N, O, S, Si and Fe. We have collected the most recent data on these galaxies and compared with our model 
results. We have also collected data for Damped-Lyman a-systems. 

Results. Our results show that in order to reproduce all the properties of these galaxies, including the spread in the 
chemical abundances, the star formation should have proceeded in bursts and the number of bursts should be not 
larger than 10 in each galaxy, and that metal-enhanced galactic winds are required. A metal-enhanced wind efficiency 
increasing with galactic mass can by itself reproduce the observed mass-metallicity relation although also an increasing 
efficiency of star formation and/or number and/or duration of bursts can equally well reproduce such a relation. 
Conclusions. Metal enhanced winds together with an increasing amount of star formation with galactic mass are required 
to explain most of the properties of these galaxies. Normal galactic winds, where all the gas is lost at the same rate, 
do not reproduce the features of these galaxies. On the other hand, a global increase of the amount of star formation 
(increasing efficiency and/or number of bursts and/or burst duration) with galactic mass is able by itself to reproduce 
the mass-metallicity relation even without winds, but without metal-enhanced winds is not able to explain many other 
constraints. We suggest that these galaxies should have suffered a different number of bursts varying from 2 to 10 and 
that the efficiency of metal-enhanced winds should have been not too high (X m w ~ 1). We predict for these galaxies 
present time Type la SN rates from 0.00084 and 0.0023 per century. Finally, by comparing the abundance patterns of 
Damped Lyman-a objects with our models we conclude that they are very likely the progenitors of the present day 
dwarf irregulars. 
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1. Introduction gas content and mostly young stellar populations. All these 

features indicate that they are poorly evolved objects, ei- 

Galaxy formation and evolution is one of the fundamental ther newly f orme d galaxies or evolving slowly over the 

problems in astrophysics. According to hierarchical clus- Hubble time. Especially BCDs, the least chemically evolved 

tering models, larger galactic structures build up and grow st ar-forming galaxies known in the universe (12+log(0/H) 

through the accretion of dwarf galaxies which are the first rang i ng between 7.1 and -8.4). are excellent laborato- 

structures to collapse and form stars ([White fc Frenk |l99l| ries for studying nucleosynthesis processes in a metal- 

iKauffmann et al. ||1993|). These building-block galaxies are deficient environment, in conditions si milar to those pre- 

too faint and small to be studied at high redshifts, while a vailing at the time of g alaxy formation (|Thuan et al. lll99l 

class of nearby metal-deficient dwarf galaxies offer a much | Izotov fc Thiia^lTooa) . The study of the variations of one 

better chance of understanding it (|Thuan ||2008D . chemical element relative to another in these poorly evolved 

Dwarf galaxies, defined arbitrarily as galaxies hav- star-forming galaxies is crucial for our understanding of the 

ing an absolute magnitude fainter than M B ~ -18 early chemical evolution of galaxies and for constraining 

mag, are the most numerou s (about 80% - 90%) galax- models of stellar nucleosynthesis. 

ies in the nearby universe (iMateo I Il998t iGrebel I l200lt „ T , . , . . . . . . , TT 

H? T. — I Ail Innn^h rn~- j r, u Dlrrs are dominated by scattered bright Hii regions in 

Karachentscv et al. 2004 . Ihcir space density has been ,. ... . .. . TT J , ° , . r ° . ... 

siKrtrested to he about 40 times rnVher than that of bripht the °P tlcal > whlle 111 Hl the y show a complicated fractal-like 

o Uggcb Led L(J Ut clLJUUL 4:U LllllCb IllJ^llLl LIld.Il Hldl Ui Ull^llL j? U 11 -CI j- J 1 T> * 1 TT 



galaxies jStavelev-Smith et al. I Hfjh . Late-type dwarfs, l^TllfM^ ^ dumpS - Typi ° al Hl m&SSeS 

dwarf irregular galaxies (dlrrs) and blue compact dwarf — ®' 

galaxies (BCDs), are galaxies harboring active or recent BCDs are currently undergoing an intense burst of star 

star formation activity, but have low metallicitics, large formation which gives birth to a large number (10 3 — 
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10 4 ) of massive stars in a compact region( < 1 kpc), 
which ionizes the interstellar medium, producing high- 
excitation super giant Hn regions and enriching it with 



excitation supcr giant rill regions and enriching it with 
heavy elements (jThuan et al. 1119 95^. Part of the extended 
neutral g as may be kinemat ically decoupled from the 
galaxies (Ivan Zee et al. I I1998D . The majority of BCDs 
(more than 99%) are not primordial systems, but evolved 
dwarf galaxies where starburst activit y is immersed within 
an old extended stellar host ga l axy (Kunth et al. 1 119881 : 
iPapaderos et al. I Il996t iThuan I 120081 ). However, recent 
work lends strong observational support to the idea that 
some among the most metal-deficient star-forming galax- 
ies known in the local universe have formed most of 
their stellar mass within the last 1 Gyr, hence they qual- 



ify as young galaxy candidates (IPapaderos et aL 
Izotov fc Thuanll2004btlPustilnik et al. 112004 1 Aloisi et alTI 



i qual- 
2001 



20071). 



ISearle et al. I (|l973l) concluded that extremely blue 
galaxies should have undergone intense bursts of star for- 
mation (SF) separated by long quiescent periods (bursting 
SF). Recent detections of old underlying stellar populations 
in most BCDs seem to corroborate their suggestion and re- 
veal at least another burst of SF besides the present one, 
even in the case of the most metal-poor BCDs known, I Zw 
18 (lOstlinll2000l ) and SBS 0335-052W dLipovetskv et al. I 



Aside from the bursting SF mode, gasping ( Tosi et al 



1991 ) or mild continuous (jCarigi et al. I 119991 : iLegrand 



20001: iLegrand et al. Il200(il) SF regimes have been proposed 



for dirrs and BCDs. The gasping scenario, in which the 
interburst periods are significantly shorter than the active 
phase s, is probably the more realis tic picture for many of 
them (jSchulte-Ladbeck et al. 1120011) . 

It is very likely that dirrs and BCDs have suffered 
galactic winds. In the past years, theorists have argued 
that winds carry heavy elements out of galaxies, and 
that they remove a larger fraction of the metals in 
lower mass galaxies 
iDe Young fc Gallagher 



dLarson 
1990t 



1971 iDekel fc Silk I [l 986: 



Mac Low fc Ferrara 1 119991 



The observational evidence of outflows from dwarf galaxies 
has grown rapidly in time ( e . g . , iMeurer et al. 1119921 : iMartin I 
Il996t iBomans et al. 1119971) . Only recently, however, has it 
become possible to directly measure the metal content of 
galactic winds and confirm that w inds are indeed metal 
enhanced (e.g.. IMartin et al. Il2002l ). 

In the past years many models for the chemical 
evolution of these galaxies appeared and tried to ex- 
plain the intrinsic spread observed in their proper- 



ties dMatteucci & Chiosi 1983 


; Matteucci & Tosi 1 19851 


Pilvugin 1993: iMarconi et al. 


19941: iBradamante et al. 1 


19981: Henrv et al. 1 [20001: Lanfranchi & Matteucci 120031: 


Romano et al. 1 20061) . Most of these papers suggested that 



the spread can be reproduced by varying the efficiency 
of star formation or galactic wind from galaxy to galaxy 
or by assuming that there is self-pollutio n in the H11 re - 
gions where the abundances are measured (jPilvugin Ill993l ) . 
Metal-en hanced winds with d i fferent prescriptions were 
studied (IMarconi et al. 11994 IBradamante et al. I Il998t 
IRecchi et al. Il200ll 12004: iRomano et al. 1120061) . " 

Besides the chemical abundances, the photometric 
and spectral properties are also taken into a c count 
in some theoretica l work s (e.g.. IVazauez et al. I 2003 : 



Stasinska fc Izotov 20031: iMartm-Manjpn et al. 



2008, 



ferent codes of chemical evolution, evolutionary popula- 
tion synthesis and photoionization, and concluded that the 
closed box models with an attenuated bursting SF and a 
initial star formation efficiency (SFE) e = 0.1 ~ 0.3 can re- 
produce the observed abundances, diagnostic diagrams and 
equivalent width-colour relations of local H 11 galaxies. 

In this work, we present a new series of chemical evolu- 
tion models for di rrs and BCDs. The model s are based on 
the original one of IBradamante et al. I (I1998D . but consider 



a larger number of chemical species and updated stellar 
yields. We have tested both the bursting and the continuous 
regime of star formation By comparing our model results 
with the most recent data of dirrs and Damped Lyman-a 
systems (DLAs), we aim at understanding the importance 
of galactic winds, the history of star formation and the ori- 
gin of the mass-metallicity relation in dwarf irregulars. A 
comparison between the properties of local dirrs and those 
of high redshift DLAs will allow us to understand the na- 
ture of DLAs and establish whether they can be considered 
as the progenitors of local dirrs and BCDs. Moreover, by 
studying in detail the abundance patterns such as [X/Fe] 
versus [Fc/H] will allow us to understand if these dwarf 
galaxies can be the building blocks of more massive galax- 
ies, as suggested by the hierarchical clustering scenario of 
galaxy formation. 

This paper is organized as follows. In Section 2, we 
present the observational constraints. In Section 3, the 
adopted chemical evolution models are described. Our 
model results are presented in details in Section 4 and dis- 
cussed in Section 5. 

2. Observational Properties 

The mctallicity, defined as the fraction of elements other 
than hydrogen and helium by mass, is an important indi- 
cator of the formation and evolutionary stage of a galaxy, 
and usually correlates with macroscopic properties of late- 
type galaxies, e.g. luminosity, mass, gas fraction, rota- 
tion speed, mo r pholo g ical type, etc. (e .g., Garnett _ 20021: 



or pholo gica 

Pilvugin et all 120041: iLee et al I 120061: IVaduvescu et al. 1 
2007D . Except hydrogen and helium, oxygen is the most 
abundant element in the universe and easy to be measured 
in H11 regions because of its bright emission lines. In prac- 
tice, the oxygen abundance is usually used to represent the 
mctallicity of the galaxy. 



2.1. Luminosity-Metallicity (L — Z) relation and 
Mass-Metallicity (Al - Z) relation 

The strong correlation between the metallicity Z and the lu- 
minosity L of a galaxy is a robust relationship, holding over 
10 mag in galax y optical luminosity and a factor of 100 in 
meta lli city (e.g.. [Garnett fc Shields"! 1987 Zaritskv et al. I 



19941: iLamareille et al. I l2004t iTremonti et al. I I2004D 



Leaueux et al. I ( 19791 ) have shown first the existence of 



the correlation between the metallicity and the mass in 
both compact a nd irr egular galaxies, then confirmed by 
ISkillman et al. I (|l989l ) who found that more luminous 
(or more massive) galaxies are more metal rich. The 
correlation is also found in spirals and elliptical galax- 



ies fe.g..|Garnett fc Shields 
Zaritskv et al. 1119941: Ivan Zee et al 



1987; Brodie fc Huchra 



1991 



. .19971: ITremonti et all 
2004l ; lLee et al. ll2006lTvaduvescu et al. 1120071) . Since the lu 



20091 ). Martin-Manion et al. I ( 20081 . l2009h combined dif- minosity of a galaxy closely relates to its stellar mass, the 
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best fit of Lee et al. (2006) 
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log M. (M ) 

Fig. 1. Oxygen abundance vs. stellar mass for nearby dlrrs 
and BCDs. Blue filled a nd m agenta open circles are data 
from IVaduvescu et al. I (|2007|) . and represent BCDs and 
dlrrs res pectively; magent a open triangles are dlrrs ob- 
served bv lLee et al. I (|2006f ). and the dashed line shows the 
best fit of their data. The solid line is the best linear fit of 
all data. 



L — Z relation should represent also the mass-mctallicity 
relation. But how the luminosity is representative of stel- 
lar mass it depends on the frequency band one investi- 



cal wavelengths (e.g., Leaueux et al. 1979: Skillman et al. 


19891 1997: Pilvuginl 


20011: Garnett 1 120021: iLee et al. 


2003allbl: iPilvugin et al. 


20041: van Zee & Havnes 1 120061: 


Ekta & Chensralur 1 l201f 


h. The optical luminosity could 



be affected by the current star formation process, there- 
fore more and more efforts were put into the determi- 
nation of the near-infrared L — Z (and hence M* — Z) 
relation where the dominant emission arises from the 
older stellar populati ons 
20041 " 



J.C. Lee et al. 



■g- : 



Perez-Gonzalez et al. I 2003 ; 
Salzer et"al 1 120051: ILee et al I 2006 : 
Mendes de Oliveira et'a l. I 120061: iRosenberg et al. I booa 



Vaduvescu et al. I 120071 : ISaviane et al. I 120081 ). Lee et al. 1 
(|2006l) considered 27 nearby star-forming dwarf irregu- 
lar galaxies whose masses spread over 3 dex, and exam- 
ined the — Z relation at 4.5 iim (Spitzer)(see Fig. [TJ. 
IVaduvescu et al. I (|2005ll2006ll2007h studied the properties 
of both dlrrs and BCDs, and obtained the M» — Z rela- 
tion by assuming M*/Lk = 0.8 Mq/Lkq- They concluded 
that, for both dlrrs and BCDs, metallicity correlates with 
stellar mass, gas mass, and baryonic mass, in the sense that 
more massive systems are more metal-rich. 



2.2. Gas fraction-metallicity (fj, — Z) relation 

A more useful relation for chemical evolution studies is the 
metallicity-gas fraction (n — Z) relation, because it provides 
the information about how the gas convert into stars and 
metals, and also about gas flows (infall and/or outflow), 



O 



C\2 




2007) 
2007) 



• BCG (Vaduvescu et al. 
O dlrr (Vaduvescu et al. 
A dlrr (Lee et al. 2006) 

* Irr (Pilyugin et al. 2004) 
■ Irr (Garnett 2002) 

closed box model (y z -0.007) 



-0.5 -1 -1.5 

log(M gas /M tot ) 

Fig. 2. Oxygen abundance vs. gas fraction relation for 
nearby dlrrs and BCDs. Blue filled and ma genta open cir- 
cles are data from IVaduvescu et al. I (|2007| ). and represent 
BCDs and dlrrs r espectively ; mage nta open triangles are 
dlrrs observed by ILee et ahl (|2006l); cyan filled tr i angles 
and squares a r e the Irrs observed bv lPilvugin et al. 1 (|2004[) 
and lGarnett I (|2002h respectively. The solid line is the pre- 
diction of the closed box model (Z = y^ln^ -1 ) with an 
effective yield yz = 0.007 correspo nding to the Salpeter 
IMF and lWooslev fc Weaver I ([l995) nucleosynthesis. 



as first suggested by iMatteucci fc Chiosi I (|1983[ ) in their 
chemical evolution models of dwarf irregular galaxies. 

As former works have shown, if the system does not have 
gas flows, the metallicity evolution predicted by the closed 
box model is a simple function of the gas fraction \x and 
true yield y unde r the instantaneous recycling assumptions , 
Z = yzHpi- 1 ) dSchmidt I [19631 : iSearle fc Sargent lll97^l . 
By comparing the "observed" effective yield, yz,eff = 
Z bs/ ln(^t~ bs ) with the true measured yield yz-, one can 
understand whether the system evolved as a closed box or 
if infall and/or outflow have been important. In fact, both 
infall and outflow have the effect of decreasing the effective 
yield; 

ILee et al. I (|2003allbT ) showed that the oxygen abundance 
is tightly correlated with the gas fraction in dwarf irregu- 
lar galaxies through optical observations, and they argued 
that these dlrrs have evolved in r elative isolation, with- 
out inflow or outflow of gas. Lately. ILee et al.l (|2006h mea- 
sured the 4.5 /im luminosities for 27 nearby dlrrs with 
the Spitzer Infrared Array Camera, and showed the re- 
lation between oxygen abundances and gas-to-stellar mass 
ratio. Their results suggest reduced yields and/or signifi- 
cant outflow rate s, which have been also indicated by prev i- 
ous authors (e.g. JGarnett Il200"2l Ivan Zee fc Havnes 1120061 ). 
Using NIR photometry , the most important discovery of 
IVaduvescu et al. I (|2007l ) is that the ii — Z relation for BC Ds 
follows that of dlrrs, and agrees with ILee et all (|2003aU bT) 
that the evolution of field dlrrs has not been noticeably 
influenced by gas flows. However, they also point out that 
several dlrrs and at least one BCD do show HI deficiencies 
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Fig. 3. Abundance ra tio of different elements. Blue open circles arc BCDs, from llzotov fe Thuan I (119991 ) and 
iPapaderos et al. I (|2006f ); green open squares arc DLAs, the data are shown in Table [TJ 



in dense environments, indicating that the gas may be re- 
moved by external processes. Using the effective yield as the 
only criterion for gas flow is too simplistic and the conclu- 
sion might not be reliable. Convincing conclusions should 
be reached through detailed modeling work. 

In Fig. [2J we show the observed oxygen abundance-gas 
fraction relation of nearby dlrrs and BCDs. ^From this 
figure we can see, the BCDs have larger abundance and 
gas fraction range, while most dlrrs have higher gas frac- 
tion compared with BCDs and Irrs, implying their poorly 
evolved stages. 

2.3. Abundance ratios 

How the abundances of chemical elements change relative to 
one another is a crucial clue for understanding the chemical 
evolution of galaxies and stellar nucleosynthesis. 

Hn regions are ionized by newly born massive stars, 
hence showing the metallicity of the ISM at the present 
time. Therefore, metallicities in dlrrs and BCDs are usu- 
ally derived from the ionized gas in Hn regions through 



their strong narrow emission lines, llzotov fc Thuan I (|l999t) 
presented high-quality ground-based spectroscopic obser- 
vations of 54 supcrgiant Hn regions in 50 low-mctallicity 
BCDs with oxygen abundances 12+log(0/H) between 7.1 
and 8.3, and determined abundances for the elements N, O, 
Ne, S, Ar, Fe, and also C and Si in a subsample of 7 BCDs. 
IPapaderos et al. I (|2006t ) presented spectroscopic and pho- 
tometric studies of nearby BCDs in the 2dFGRS (Two- 
Degree Field Galaxy Redshift Survey) , and measured their 
Ne/O, Fe/O and Ar/O ratios. Both of these works do not 
consider the dust de pletio n correction. We show the data 
of llzotov fe Thuan I (|l999t) and IPapaderos et al. I (|2006t) in 

Fig. m 

DLA absorption systems, found in the spectra of high- 
redshift QSOs, arc neutral clouds with large Hi column den- 
sities, JV(HI) > 2 ■ 10 20 cur 2 . They arc likely to be proto- 
galactic clumps embedded in dark matter halos and may 
provide the important information on the early chemical 
evolution of galaxies. With high resolution spectroscopy 
of QSO absorption lines, elemental abundances can be 
measured up to redshift z w 5. The chemical abun- 
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dances of DLA systems give us complementary observa- 
tional constraints on the formation and evolution of galax- 
ies. Abundance measurements in DLA systems relevant for 
the present work are listed in Table [1] and plotted in Fig. [3] 
(green squares). In comparing DLA abundances with model 
predictions care must be taken for dust depletion effects. 
Luckily, these effects are expected to be negligible for most 
of the elements used in the present investigation, such as 
C, N, O and S: in fact these elements show litt le values 
of de pletion, if any, in nearby interstellar clouds (jjenkins I 
|2009() and are expected to be even less depleted in DLA 
systems. On the other hand, we expect some depletion ef- 
fects for Fe and, to a lesser extent, for Si. Estimates of Fc 
depletion i n DLAs based on the comparison with Zn mea- 
surements (|Vladilo I I2004T ) are available only for a few sys- 
tems of Table [T] These results indicate that Fe tend to be 
underestimated when the level of metallicity is relatively 
high. This explain the few cases with largest deviations 
from BCD measurements and from the model predictions 
shown in Figs. H El [HI UM CCD HI and \T7\ 

2.4. Primordial helium abundance Y p and AY/AZ 

The determination of primordial helium abundance, Y p , is 
important for the study of cosmology and the evolution 
of galaxies, because an accurate initial Y is required to 
test Big Bang nucleosynthesis and build chemical evolution 
models. 

One way of estimating Y p is by extrapolating the ob- 
served helium-metallicity (Y — Z) relation to Z = by 
assuming the slope AY/AZ to be constant. More re- 
cently, it has been common practice to use AY / AO since 
the oxygen abundance is easier to determine and can 
represent the metals. To obtain an accurate Y p value, 
a reliable determinatio n of AY / AO for oxygen-poor ob- 
jects i s needed (e.g., Ilzotov et al. I 19991 A. Peimbert I 



20031: iLuridiana et al. I 120031: Ilzotov fc Thuan I l2004al; 



M. Peimbert Il2007f ). Ilzotov fc Thuan I (|2004al) derived the 
primordial helium Y p = 0.2429 ± 0.0009 and the slope 
AY / AO = 4.3 ± 0.7 from observations of 82 Hn regions. 
For a restricted sample (7 Hn regions), they obtained 
F p = 0.2421 ± 0.0 021 and AY / AO = 5.7 ± 1.8. Later, 
Ilzotov et all (|2006fl derived Y p = 0.2463 ± 0.0030 from the 
emission of the whole Hn region of the extremely metal- 
deficient blue comp act dwarf galaxy SBS 0335 — 052E. 
IM. Peimbertl (|2007h has adopted AY / AO = 3.3 ± 0.7 
from theoretical and observational results, and derived 
Y p = 0.2474 ± 0.0029. These va lues are in exc e llent agree- 
ment with the value derived by ISoereel et al. I (|2007h from 
the WMAP results, Y p = 0.2482 ± 0.0004. 

In Fi g. [H we replot th e helium -oxygen abundance re- 
lation of Ilzotov fc Thuan I (|2004ah by* using their data in 
Table 5. The linear regression is the one derived from the 
whole sample Y = 0.2429 + 43 * (O/H). 



3. Model Prescriptions 

In this work, we used an up dated version of the chem ical 
evolution model developed bv lBradamante et al. I (|1998t) to 
study the formation and evolution of late-type dwarf galax- 
ies, dlrrs and BCDs. 

The general picture is the following: our model is one- 
zone and assumes the galaxy built up by continuous infall 
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from Ilzotov et al. (119971): magenta open pentagrams , from 
other data (llzotoy fc Thuan I Il998at Hzptoy et al. I 119991: 



Thuan et al 1 119991: Ilzotov et al. I l2001al lbl: iGuseva et al. I 
2001L I2003allbh : the solid line is the maximum likelihood 
linear regression of all data, Y = 0.2429 + 43 * (O/H) 



of primordial gas (X = 0.7571, Y p = 0.2429, Z = 0). Stars 
form and then contaminate the interstellar medium (ISM) 
with their newly produced elements which mix with the 
ISM instantaneously and completely. Stellar lifetimes are 
taken into account in detail, i.e. the instantaneous recycling 
approximation (IRA) is relaxed. The energy released by su- 
pernovae (SNe) and stellar winds is partially deposited in 
the ISM, and galactic winds develop when the thermal en- 
ergy of the gas exceeds its binding energy. The wind expels 
metals from the galaxy, hence it has a significant influence 
on the chemical enrichment of the galaxy. 

The time evolution of the fractional mass of the element 
i in the gas, Gi, is described by the equations: 

Gi(t) = -ip(t)Xi(t) + Ri(t) + G iMf {t) - G hOUt (t), (1) 

where Gi(t) = M g {t)Xi(t) /M^tc) is the gas mass in the 
form of an element i normalized to the total baryonic mass 
Ml at the present day to = 13 Gyr; M g (t) is the gas mass 
at time t and Xi(t) represents the mass fraction of ele- 
ment i in the gas, i.e., abundance by mass. The quantity 
G(t) = Mg(t)/ML(tG) represents the total fractional mass 
of gas and Xi(t) can be expressed by Gi(t)/G(t). The four 
items on the right hand side of equation (JTJ) show the mass 
change of the element i caused by the formation of new 
stars ip(t)Xi(t), the material returned through stellar winds 
or SN explosion Ri(t), the infall of primordial gas Gi^ n f(t), 
and the outflow Gi tOU t(t) respectively. 

The star formation rate (SFR) ip(t) in this work is sim- 
ply assumed as: 

m = eG(t), (2) 
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Fig. 5. The evolutionary tracks of [O/Fc] vs. [Fc/H] as predicted by models without outflow, a) models with different 
number of bursts; b) models with different occurrence times of bursts; c) models with different durations of bursts; d) 
models with different star formation efficiency. All these models assume the same total infall mass Mj„f = 10 9 M Q . The 
cyan open circles are BCDs and green open squares are DLAs, same as in Fig. [3] 



where e is the star formation efficiency and is in units of 
Gyr , being one of the free parameters in our work. 

The rate of gas infall is assumed to be exponentially 
decreasing with time: 

f mf (t)=Ae- t /\ (3) 

where A is the normalization constant which is constrained 
by the boundary condition J Q G Ae~'/ T = 1, and r is the 
infall timescale. So we can easily obtain the accretion rate 
of an element i through the formula 

GiSnf(t) — Xi^ n ffi n f — — ■ '— ^_ tc ^ T ^ , (4) 

inf — 0(* 7^ H, He) if primordial gas is assumed. 
In our model, the galactic wind develops when the ther- 
mal energy of the gas Eg h (t) exceeds its binding energy 

E b g (ty. 

Ef{t) > E b g (t). (5) 



The thermal energy of the gas is produced by SN explosions 
(both Type II and Type la) and stellar winds: 

Ef{t) = Ef NII (t) + Ef NIa (t) + Et h w (t). (6) 

However, not all the energy produced in the above men- 
tioned events is stor ed into the ISM, since a f raction of it 
is lost by cooling. In iBradamante et al. I (|l998T )'s work, the 
efficiencies of energy transfer from SN and stellar winds 
into the ISM are the same, wsnii = VSNia = Vsw = 0.03 
(see their work for more details). However, m ore recently, 
iRecchi et al~l (|200lD and IRecchi et al~l (|2002D have shown 
that since SN la explosions occur in a hotter and more 
rarefied medium, their energy can be more efficiently ther- 
malized into the ISM and, consequently, their efficiency of 
energy transfer is higher. Therefore, in this work we as- 
sume that the efficiencies of energy transfer are t]snii = 
0m,v S Nia = 0.8 and n sw = 0.03 for SN II, SN la and 
stellar winds respectively. 
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Fig. 6. Present oxygen abundance vs. gas fraction {upper 
panel) and Y vs. oxygen abundance (lower panel) as pre- 
dicted by models without wind. All these models assume 
the same total infall mass Mi n f — 10 Mq. Big black open 
circles: models with different numbers of bursts; big black 
open triangles: models with different durations of bursts; big 
black open squares: models with different occurrence times 
of bursts; big black open pentacles: models with different 
star formation efficiency. Black solid, red dotted and blue- 
dash lines in upper panel are the evolutionary tracks of 
model M9d4, M9tau and M9IMF. The black-solid line in 
lower panel is the best fit of all the model points. The ob- 
servational data are the same as in Fig. [2] and Fig. [H cyan 
open circles and cyan open triangles represent BCDs and 
dlrrs respectively. The cyan solid line is the best fit of the 
observational data. 



To c ompute Eg(t), the b inding energy of gas, we also 
followed iBradamante et al~l (|1998l) and assumed that each 
galaxy has a dark matter halo. The binding energy of gas 
is described as: 



E B gas{t) = W L (t) + W LD (t) 



with: 



W L {t) 



0.5 G 



M gas {t)M L {t) 



I'L 



(7) 



(8) 



which is the potential well due to the luminous matter and 
with: 

M gas (t)M dark 



W LD (t) = — Guild - 



(9) 



which represents the potential well due to the interac- 
tion between dark and luminous matter, where wld ~ 
j^S(l + 1.37S), with S = rj^jrY), being the ratio between 
the galaxy effective ra dius fry,) and the ra dius of the dark 
m atter core (rp) (see iBertin et al.1ll992[ ).We assumed as 
in IBradamante et al~1 ( 19981 ) that the dark matter halo is 
10 times more massive than the luminous matter and that 
5 = 0.3. 

The rate of gas loss via galactic wind for each element 
is assumed to be simply proportional to the amount of gas 
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Fig. 7. Present day oxygen abundance vs. M* as predicted 
by models without wind. The upper (middle/lower) panel 
shows the model results with different numbers of bursts 
(durations/SFE), along the blue- dashed lines the number 
of bursts (durations/SFE) increases from 2 to 7 (0.02 to 
0.3/ 0.1 to 2), and along the black-solid line the total infall 
mass Mi n f increases from 10 8 to 1O 1O M0. The red-thick- 
solid line connects the models which are consistent with the 
observational data, larger n (d/e) for more massive galaxy. 
The observational data are the same as in Fig.[TJ cyan open 
circles for BCDs and cyan open triangles for dlrrs. 



present at the time t: 

Gi, ou t(t) 



Wi\G(t)Xi tOUt {t), 



(10) 



where Xi tOUt (t), the abundance of the element i in the wind, 
is assumed to be same with -Xj(t), the abundance in the 
ISM; A describes the efficiency of the galactic wind and has 
the same units as e (Gyr -1 ); u>i is the efficiency weight of 
each element, hence WiX is the effective wind efficiency of 
the element i. X and Wi are the other two free parameters 
in our model. In this work, we have studied two kinds of 
wind, the normal wind and the metal-enhanced wind. In the 
case of the normal wind, all elements are lost in the same 
way, i.e., ittj = 1 for all elements, and we use A^, to denote 
the wind efficiency in this case; however, wind provoked 
by SN explosion could carry out more metals than H and 



enhanced" wind dMac Low & Ferrara 


19991: iRecchi et al. 1 


l200lt iFuiita et al. 20031: Recchi et al. 


20081). and we use 



A 



v to denote the wind efficiency in this case. 
The initial mass function (IMF) is usually assumed to 
be constant both in space and time in different galaxies, 
and can be expr essed as a powe r law of stellar mass as 
suggested first bv ISalpeter I (|l955h : 



-(1+x) 



(11) 



where x = 1.35 for Salpeter IMF and </>o is the nor- 
malization constant which can be obtained by satisfying 
J m<f>(m)dm — 1 in the mass range 0.1 — 100 Mq. However, 
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Table 2. Parameters of models without outflow 
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Fig. 8. The evolutionary track as predicted by models with normal wind. The left two panels are log(C/0) vs. 
12+log(0/H) and [O/Fc] vs. [Fe/H]; the right two panels are oxygen abundance vs. gas fraction and Y vs. (O/H). 
All the 5 models have same bursts sequence (t = 1 Gyr, t = 1/3/5/7/9/11/13 Gyr, and d = 0.1 Gyr for each burst), 
but different wind efficiencies, green-solid lines for A TO = 0, black- dash- dot lines for A TO = 0.2, red-long-dash lines for 
A m = 0.5, and blue- short- dash lines for A m = 1. A model with long infall timescale (r = 10 Gyr, A m = 1) is shown in 
magenta dotted line. All these models assume the same total infall mass Mi n f = 10 9 M Q . The observational data are the 
same as in Fig. [21 Fig. [3] and Fig. SI BCDs, dlrrs, and DLAs are plotted in cyan open circles, cyan open triangles and 
green open squares here. 
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we tested also the lScalol (jl986ft IMF: 



4>scaioS&(rn) cx 



-2.35 



-2.7 



(0.1 s$ m < 2) 
(2^m^ 100). 



(12) 



which is a two-slope IMF and is steeper at the high mass 
end. 

Stellar yields of different elements are important ingre- 
dients of chemica l evolution studies. In this w ork, we adopt 
stellar yie lds of Wooslev fc Weaver I (119951) fo r massive 
stars and Ivan den Hoek fe Groeneweeen I ( 1997T ) for low- 
and intermediate-mass stars. Both of them are metallicity- 
dependent. 



4. Model Results 

In order to understand the observed global proper- 
ties and abundance patterns of late-type dwarf galax- 
ies, we have calculated several models. The typical 
galaxy is assumed to be forming by continuous in- 
fall of primordial gas and bursting star f ormation, as 



19731 Matteucci & Chiosi 1 19831 Marconi et al. I 


1994; 


Bradamante et al. 1 


1998; Lanfranchi & Matteucci 1 


2003; 


Romano et al. 2006 


) . We also checked the case of continu- 



ous star formation for these galaxies (see Sect. 4.5). 



4.1. Model without outflow 

Here we examine the models without outflow. Different 
numbers n, durations d, times of the occurrence of bursts t, 
and different star formation efficiencies e have been tested. 
The parameters adopted in the models are listed in Table [5] 
All the models assume short infall timcscalcs (r = 1 Gyr) 
except model M9tau, and the galactic lifetime is taken to 
be 13 Gyr for all the models. From the second to the forth 
columns there are the model names, classified by different 
total infalling mass, from 1O 8 M to 10 10 M Q ; the fifth col- 
umn shows the SFE; the sixth column the number of the 
bursts; the seventh column the middle time of each burst; 
the eighth column the duration of each burst, where " 0.1*3 
" means the duration of three bursts are the same, namely 
0.1 Gyr. 

In Fig. [SI we show the evolutionary tracks of [O/Fe] 
vs. [Fe/H] as predicted by models with a different num- 
ber of bursts(panel a), different times for the occurrence of 
the bursts (panel b), different burst durations (panel c), and 
also different SFEs (panel d). The most distinctive feature 
of the bursting star formation scenario is the "saw-tooth" 
behaviour of the tracks, which is caused by the different 
origins of en-elements and iron-peak elements. Oxygen is 
mainly synthesized by massive stars, therefore its abun- 
dance increases only during the bursting time, whereas iron 
is mainly synthesized by the Type la supcrnovae and its 
abundance still increases after the burst is over owing to 
the time delay of SN la explosions. 

As we can see from Fig. [5j different durations of bursts 
and different SFEs among galaxies could be the explanation 
of the scatter in the data. In panel c, as a comparison, 
we also plot the model results for a low continuous star 
formation process (model M9d4). Clearly in this case the 
saw-tooth behaviour disappears. 
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Fig. 9. Present day oxygen abundance vs. M* as predicted 
by models with normal wind. All the models contain 7 
bursts and the same burst sequence (e = 0.2, n = 7, t = 
1/3/5/7/9/11/13, d= 0.1 for each burst), and the total in- 
fall mass varies from Mj„/ = 10 s to 10 10 Af Q . The upper 
panel shows results with short infall timescale (r = 1 Gyr) , 
whereas the lower one shows results with long timescale 
(r = 10 Gyr). The predicted M — Z relation with dif- 
ferent strength of normal wind are shown in black lines: 
black solid, dash-dot, long-dash, short-dash, and dotted lines 
are for A TO = 0, 0.2, 0.5, 1, 3 respectively. Models with same 
Mi n f are connected by blue lines. An increasing wind ef- 
ficiency to less massive galaxies is shown in red-solid line. 
Observational data are shown in cyan open circles (BCDs) 
and cyan open triangles (dlrrs), the cyan solid line is the 
best fit of all data, same as in Fig. [TJ 



If we compare the present oxygen abundance vs. the gas 
fraction as predicted by the above models with the obser- 
vations (Fig. [6j upper panel), we can see that no matter 
how the star formation history changes (different n, t, d, e), 
the model results always stay along the sam e curve. As 
first suggested bv iMatteucci fc Chiosi I (|l983l ). in order to 
explain the spread in this diagram one should necessarily 
claim a variation of the IMF, or of the wind rate, or of the 
infall rate. Therefore, we examined the other IMF (jScalo I 
119861 ) which is steeper than Salpeter IMF at the massive 
end. However, the model results in Fig. [5] still cannot ex- 
plain the observed lowest oxygen abundance at the same /x 
or the lowest /j, at the same oxygen abundance. Therefore, 
unless one assumes unrealistically steeper IMFs, the ob- 
served fi — Z strongly implies that there should have other 
mechanisms operating in the galaxy which can reduce the 
O abundance or the gas fraction. 

In the lower panel of Fig. [51 it is shown the Y— (O/H) 
relation as predicted by our model which is consistent with 
the observations, although a little flatter than the observed 
best fit. The best fit to these model results is Y = 0.2428 + 
12.07(O/H). 
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Fig. 10. The evolutionary track as predicted by models with various amounts of metal enhancements. The left two panels 
are log(C/0) vs. 12+log(0/H) and [O/Fe] vs. [Fe/H]; the right two panels are oxygen abundance vs. gas fraction and Y 
vs. (O/H). All the 3 models have same bursts sequence (t = 1/3/5/7/9/11/13 Gyr, and d = 0.1 Gyr for each burst), but 
different degrees of metal enhancement. The black-solid lines, red-dash lines, blue-dot lines represent u>n,He = 0.1,0.5, 1 
respectively. All these models assume the same total infall mass Mj n j = 10 9 M Q . The observational data are the same as 
in Fig. [51 Fig. [3] and Fig. 01 BCDs, dlrrs, and DLAs are plotted in cyan open circles, cyan open triangles and green open 
squares here. 



Since we assume a linear correlation between star forma- 
tion rate and gas, the models with different total infalling 
masses are self-similar. Therefore, when the abundance ra- 
tios and gas fraction are examined (Figs.[S]and[H]), the three 
scries of models (M8, M9, and M10) are overlapping. 



In Fig. [7]we plot the mass-metallicity relations predicted 
by models without wind. We run models for three differ- 
ent infall masses (10 s , 10 9 , 10 10 M o ). The effects of dif- 
ferent numbers of bursts (n — 1,3,7), different durations 
(d = 0.03, 0.1, 0.3) and different SFEs (e = 0.2, 0.5, 1.0, 2.0) 
as functions of galactic mass are shown. It is evident from 
Fig. [7] that our models can very well reproduce the M-Z 
relation even without galactic wind but just assuming an 
increase of the number, or duration of bursts, or the effi- 
ciency of SF. 



4.2. Model with normal wind 

If the galactic wind has the same chemical composition as 
the well-mixed ISM, i.e. uii = 1 for all the elements, we call 
it "normal wind" . 

In Fig. [51 we show the the evolutionary tracks predicted 
by models with normal wind; abundance ratios of log(C/0) 
vs. 12+log(0/H) and [O/Fc] vs. [Fe/H] are on the left side, 
while n — Z and Y — Z relations on the right side. The 
models have the same total infall mass (Mi n f — 1O 9 M ) 
and same bursts sequence (t = 1/3/5/7/9/11/13 Gyr, with 
d = 0.1 Gyr for each burst), but different wind efficiencies 
(X w = 0, 0.2, 0.5, 1.0). Oxygen is produced by massive stars, 
therefore no oxygen will be ejected into the ISM after star 
formation ceases. On the other hand, elements, such as C 
and N produced by low- and intermediate-mass stars, and 
Fe mainly produced by SN la explosion, are continuously 
polluting the ISM after the star formation stops, owing to 
their long lifetime. Therefore, the decrease of the mass of 
gas (i.e., H and He) and the a-elements lost with the wind 
will result in a dramatic increasing of the abundance of 
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Fig. 11. The evolutionary track as predicted by models with various metal-enhanced wind efficiencies. The left two 
panels arc log(C/0) vs. 12+log(0/H) and [O/Fc] vs. [Fc/H]; the right two panels are oxygen abundance vs. gas fraction 
and Y vs. (O/H). All the 5 models have same bursts sequence (t = 1/3/5/7/9/11/13 Gyr, and d = 0.1 Gyr for each 
burst), but different wind efficiencies. Green-solid lines: X mw = 0; black- dash-dot lines: X mw = 0.2; red-long-dash lines: 
X m w = 1; blue- short- dash lines: X mw = 3; magenta-dot lines: X mw = 10. All these models assume the same total infall 
mass Mi n f = 10 9 M Q . The observational data are the same as in Fig. [21 Fig. [3] and Fig. [H BCDs, dlrrs, and DLAs are 
plotted in cyan open circles, cyan open triangles and green open squares here. 



the "time-delayed" elements. The stronger the wind, the 
higher the abundance of C or Fe relative to O predicted by 
the models. 

The main effect of normal winds is to decrease the gas 
fraction with a smaller effect on the O/H abundance, as we 
can see from the 12+log(0/H)-/j, relation (upper right panel 
of Fig. [5]). This means that models with normal wind can- 
not explain the whole spread in O /H observed at a given \i 
for these galaxies. There are two possible reasons for that. 
One is the same wind efficiency (i.e., WiX w ) for both oxy- 
gen and hydrogen in the normal wind, so that both O and 
H decrease at the same time. The other one is the short 
infall time scale assumed (r = 1 Gyr). In this case, no pri- 
mordial gas falls into the galaxies to dilute the ISM at late 
evolutionary times. Therefore, we also developed a model 
with long infall time scale (t = 10 Gyr, magenta dotted 
lines in Fig. [5]). It is clear that in this model the metal- 
licity decreases in the interburst time. Actually, the infall 
of primordial gas (i.e., H and He) results in a lower mass 
loss rate of H and He than metals, similar to the metal- 



enhanced wind case which will be further discussed in the 
next section. A very strong normal wind (e.g. X w > 0.5) 
seems unlikely in late-type dwarf galaxies since it would 
lose a large amount of gas, and hence it would predict a 
too low gas fraction, as it is evident in Fig. [5J In the lower 
right panel of Fig. [8j the predicted Y vs. (O/H) relation is 
shown and it is consistent with the observational data at 
the low metallicity, because the wind does not develop yet 
when the galaxy is still very metal poor. However, after the 
wind, an increase of the helium abundance as well as of the 
abundances of elements produced on long timescale occurs, 
especially in the case of a strong wind which produces a 
very small final gas fraction. 

In the last section, we have pointed out that the ob- 
served M — Z relation could indicate more star bursts 
or longer duration of each burst or higher SFE in more 
massive galaxies, if no outflow takes place. Now we exam- 
ine the possibility of the normal wind being the explana- 
tion of the M — Z relation, as suggested by many previ- 
ous authors. In Fig. [HI we take the models with 7 bursts 
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for example (i.e., Model M8n4, M9n4, M10n4 in no- wind 
case) but different strengthes of normal wind are introduced 
(A w = 0, 0.2, 0.5, 1, 3). It is evident from the upper panel of 
Fig. IHl that by varying only the efficiency of a normal wind 
one cannot reproduce the M — Z relation, unless other pa- 
rameters, such as the efficiency of SF or the number of 
bursts, are assumed to vary as functions of the galactic 
mass. When a long infall timescale is adopted (r = 10 Gyr, 
lower panel of Fig. [H]), less stars are formed in each model 
due to the slow gas accretion process, and the wind could 
not be induced in the high mass systems (M in f w 10 10 M Q ). 
However in the low mass galaxy (Mi n f = 10 s M®), where 
the wind could develop, the newly infalling gas dilutes the 
ISM effectively. Therefore, the general trend of the M — Z 
relation could be reproduced by combining the normal wind 
with a slow accretion process. 

In summary, the normal wind can strongly reduce the 
gas fraction but it cannot reduce sensibly the O/H. To ex- 
plain the spread observed in O /H at the same [i we should 
invoke other mechanisms, such as a continuous supplement 
of primordial gas or a mctal-cnhanccd wind (see Sect. 4.3), 
both of which imply a lower mass loss rate for H and He 
relative to metals. 

4.3. Model with metal-enhanced wind 

The galactic wind, mainly induced by SN explosion, could 
blow preferentially the metal-enriched gas out of the galaxy, 
which means metals are lost more efficiently than the gas (H 
and He). We define the wind "metal-enhanced" when the 
abundances of metals it carries out are higher than in the 
ISM. Metal-enhanced winds ha ve been already suggested 
by several dynamical wor ks (e.g. lMac Low fc Ferrara Ill999t 
iRecchi et al. |[200ll . l2002f ). In our models we simply assume 
a higher wind efficiency weight Wi for heavy elements than 
H and He. In particular, we adopt w; = l(i ^ H,He), 
W H ,Hc < 1. 

The wind models with various amounts of metal en- 
hancements are shown in Fig. 1101 All the models have same 
input parameters (e = 0.5, 7 bursts and the duration is 0.1 
Gyr for each one) except for wn.He- The model experienc- 
ing a highly enriched wind (wh,Hc = 0.1) loses very little 
gas. We show also models with mild metal-enhanced wind 
(toh,Hc = 0.5) and normal wind (u>H,He = wo = !)■ The 
evolutionary tracks for the abundance ratios show a loop 
if the wind is metal-enriched. When the wind starts, oxy- 
gen is lost more efficiently than hydrogen, hence the oxygen 
abundance within the galaxy decreases with time in the in- 
terburst phase, and elements such as C and Fe will show 
increasing abundances relative to oxygen owing to their de- 
layed restoration into the ISM. This trend continues until 
the new burst occurs. Because of the newly produced oxy- 
gen supplied to the ISM the O abundance increases and, as 
a consequence, the abundances of other elements relative to 
oxygen decrease. Therefore, the evolutionary track shows a 
loop. The lower the tUH,He, the more the metals lost, the 
lower the value that the O abundance reaches. 

The mctal-cnhanccd wind has also a dramatic influence 
on the (a— Z relation, as we shown in the upper right panel of 
Fig. 1101 The normal wind mainly reduces the gas fraction 
rather than the abundance, whereas the metal-enhanced 
wind is very powerful in reducing the metallicity of the 
galaxy. 
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Fig. 12. Present day oxygen abundance vs. M» as predicted 
by models with metal-enhanced wind (u>h.Hc = 0.1). All the 
models contain 7 bursts and the same star formation his- 
tory (e = 0.2, n = 7,t= 1/3/5/7/9/11/13, d = 0.1 for each 
burst), and the total infall mass varies from Mi n f = 10 8 
to 10 10 M o . The predicted M — Z relation with different 
strength of metal-enhanced wind are shown in different 
lines: black dash-dot, long-dash, short-dash, and dotted lines 
are for X mw = 0, 0.2, 1, 3 respective. Models with same Mj„/ 
are connected by blue lines. An increasing wind efficiency 
to less massive galaxies is shown in red-solid line which can 
fit the data very well . Observational data are shown in 
cyan open circles (BCDs) and cyan open triangles (dlrrs), 
the cyan solid line is the best fit of all data, same as in 

Fig.m 

In Fig. 111! models with different strengths (X mw = 
0,0.2,1,3,10) of highly enriched winds (wu,H e = 0-1) are 
shown. We plot the evolutionary tracks of models with same 
bursting history as in Fig. 1101 A stronger wind not only re- 
duces the abundances and increases the ratio between the 
long recycling term elements and the short ones, but also 
decreases the gas fraction dramatically. It is worth to point 
out that although the wind efficiencies X mw adopted here 
are much higher than the ones of normal wind case A^ , the 
gas is lost less effectively. To compare the results of nor- 
mal and metal-enhanced wind models, one should assume 
a larger X mw for the latter case. For example, \ mw = 10 for 
metal-enhanced winds is then multiplied by wh.Hc = 0.1, 
therefore it is comparable with the case normal wind and 
X w = 1. 

In the lower right panels of both Fig. [TU] and Fig. [TT] 

we show the Y— (O/H) relations of galaxies with different 
degrees of enriched winds u>h,Hc and different wind efficien- 
cies Xmw, but the same formation histories. By comparing 
with Fig. [51 the present-day oxygen abundances are lower 
as we expect. In this scenario, a very high helium abun- 
dance can be reached at low metallicity level, especially 
when the wind is very strong, because most of it stays inside 
the galaxy while heavy elements are lost. The present-day 
Y— (O/H) relation predicted by these models do not stay on 
a straight line in the low metellicity region, even if the unre- 



13 



J. Yin et al.: Chemical Evolution of dlrrs and BCDs 
Table 3. Parameters of the best models (u>H,He = 0.3). 





Model name 




SFE 


n 


t 1 




d 2 (Gyr) 




10 s M 


10 9 M o 


1O 1O M 


Gyr" 1 


Gyr" 1 




Gyr 


10 s M 


10 9 M Q 


10 10 M Q 


M8bl 


M9bl 


MlObl 


0.8 


0.5 


3 


3/9/13 


0.1*3 


0.3*3 


0.9*3 


M8b2 


M9b2 


M10b2 


0.8 


0.5 


5 


1/3/7/10/13 


0.1*5 


0.3*5 


0.9*5 


M8b3 


M9b3 


M10b3 


0.8 


0.5 


7 


1/3/5/7/9/11/13 


0.1*7 


0.3*6 


0.9*7 


M8b4 


M9b4 


M10b4 


0.8 


0.5 


9 


1/2.5/4/5.5/7/8.5/10/11.5/13 


0.1*9 


0.3*9 


0.9*9 




DLA1 




0.8 


0.5 


5 


1/3/7/10/13 




0.01/0.2*4 






DLA2 




0.8 


0.5 


5 


1/3/7/10/13 




0.02/0.2*4 






DLA3 




0.8 


0.5 


5 


1/3/7/10/13 




0.05/0.2*4 






DLA4 




0.8 


0.5 


5 


1/3/7/10/13 




0.1/0.2*4 





Notes. W The middle time of the burst 




12+log(0/H) [N/H] 

Fig. 13. The evolutionary track of abundance ratios of C/O, N/O, S/O, Si/O, C/Fc, N/Fc, O/Fe, and N/Si as predicted 
by our best models (Mi n f = 10 Mq ) with metal-enhanced wind (wh.Ho = 0.3). The red-solid, black- dash- dot, magenta- 
long-dash and blue-short-dash lines are the results of M9bl, M9b2, M9b3 and M9b4 in Table. El The observational data 
are same as in Fig. [31 cyan open circles for BCDs and green open squares for DLAs. 

alistic models (very strong wind X mw = 10 cases) are ruled constraints. Therefore, the observed scatter of Y— (O/H) re- 
out considering their disagreement with other observational lation may be caused by metal-enhanced winds. Based on 



14 



J. Yin et al.: Chemical Evolution of dlrrs and BCDs 



Table 4. The maximum and present-day values of the SN 
la rates by number as predicted by the best models. 



Model 
name 


number 
of burst 


rsNia (century 1 ) 


maximum 


present day 


M9bl 


3 


0.0218 


0.0016 


M9b2 


5 


0.0229 


0.0012 


M9b3 


7 


0.0229 


0.0015 


M9b4 


9 


0.0237 


0.0016 
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Fig. 14. The evolutionary tracks of 12+log(0/H) vs. gas 
fraction (upper panel) and Y— (O/H) relation (lower panel) 
as predicted by our best models (M in f = 1O 9 M ) with 
metal-enhanced wind (u>H,He = 0.3). The red-solid, black- 
dash-dot, magenta-long- dash and blue- short- dash lines are 
the results of M9bl, M9b2, M9b3 and M9b4 in Table.H The 
observational data are same as in Fig. [2] and Fig. HJ cyan 
open circles for BCDs and cyan open triangles for dlrrs. 



our model predictions, we suggest to fit the lower envelop of 
the observational data, when one derives the primordial he- 
lium Y p , because it may not be affected by the wind, hence 
the extrapolation to Z = will be more close to the real 
Y p . 

Fig. [TU is the same as Fig. |H] but for models with 
metal-enhanced wind. Unlike in the normal wind case, the 
metal-enhanced one is very effective in reducing the oxy- 
gen abundance, hence in creating the M — Z relation. 
The stronger the wind efficiency, the steeper the predicted 
AI — Z. Therefore, models with an increasing wind efficiency 
to less massive galaxies are consistent with the observa- 
tions very well. As a conclusion, the observed M — Z re- 
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Fig. 15. The SN la rate normalized to the galaxy stellar 
mass as a function of time as predicted by our best models 
(Mi n f = lO 9 -/!/©). The red-solid, black- dash- dot, magenta- 
long-dash and blue- short- dash lines are for model M9bl, 
M9b2, M9b3 and M9b4 respectively. As a comparison, the 
shaded area shows the observed range of normalized SN 
la rate i n the Irr galaxies, 0.7 7^n per century and per 
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Fig. 16. The mass-mctallicity relation as predicted by our 
best models. The total infalling masses range from Mi n f = 
10 8 to 10 10 M Q ), and more massive one prefers longer dura- 
tion of burst. The results with different numbers of burst are 
shown in different lines, red-solid, black-dash-dot, magenta- 
long-dash and blue- short- dash lines are for n = 3,5,7,9 
respectively. The observational data of BCDs (cyan open 
circles) and dlrrs (cyan open triangles) are the same as in 

Fig.rjj 



lation could be caused by metal-enhanced winds with mild 



strength (e.g., X mw < 1). 
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In summary, metal-enhanced winds should take place in 
late- type dwarf galaxies, and they play an important role 
in removing gas, especially the metals, out of the galaxies. 

4.4. Best models 

As we have shown in the last section, highly enriched or 
very strong winds will reduce the galactic oxygen abun- 
dance to a very low value during the interburst time which 
has not been confirmed from the observational point of 
view. Thus, in our best models wh,Hc = 0.3 and X mw = 0.8 
are assumed. Different numbers of bursts are examined for 
Mi n f = 1O 9 M galaxies, and the best fit models are shown 
in Fig. [T3] and Fig. [TU In these models, the same SFEs 
(e = 0.5) and burst durations (d=0.3) are assumed, and 
the details of each model are listed in Table. [3] 

We show the abundance ratios of different elements rel- 
ative to oxygen or iron as predicted by our best models for 
Mi n f = 1O 9 M in Fig.[13j The evolutionary tracks of these 
models can pass through the DLA data at early time, and 
cover the regions where most of BCDs have been observed. 

The evolutionary tracks of the fj, — Z relation (upper 
panel) and Z — Y relation (lower panel) predicted by our 
best models are shown in Fig. Q31 After the wind devel- 
ops, more are the bursts that a galaxy suffers, more are 
the oscillations in the evolutionary tracks, and these tracks 
pass through most of the data, thus explaining the observed 
spread. On the other hand, in the panel showing Y vs. 
(O/H) helium keeps increasing while O is oscillating and 
this is due to the fact that this element is produced also 
during the interburst periods and lost less effectively than 
O. This implies again that we should use the lower envelop 
of the observational data to derive the primordial abun- 
dance of helium. 

In Fig. [TSl we show the evolution of the Type la su- 
pernova rates predicted by our best models. The rates in 
this figure are normalized to the galaxy stellar mass at 
that time, i.e., expressed in number of SNe per century 
and per 1O 1O M . The peaks are always associated with 
the star formation periods, but SNe la also explode dur- 
ing the interburst times. Compared to the observed value 
quoted for the Irr galaxies, 0.77+^31 per century and per 
1O 1O M (iMannucci et al. ll2005h , our model predicts a very 
high normalized SN la rate during the first star formation 
burst, but ^1 — 2 Gyr later, the rate decays, and becomes 
comparable to the Irr galaxies in the following SF periods. 
In these best models, the maximum of SN la rate by num- 
ber varies between 0.022 and 0.024 per century and the 
present value varies between 0.0012 and 0.0016 per century 
(see Table H]). ISullivan et al~l (|2006l) have studied the rela- 
tion between SN la rate and the stellar mass of the host 
galaxy in the redshift range 0.2 — 0.75, and they found the 
SN la rate is less than 0.01 per century for star-forming 
galaxies whose stellar mass around 1O 8 M , in agreement 
with our predictions. 

In order to reproduce the mass-mctallicity relation, we 
have studied galaxies of different masses, with more massive 
galaxies preferring longer star formation bursts, d = 0.9 
Gyr for M mf = 10 10 M S and d = 0.1 Gyr for M inf = 
1O 8 M (see details in Table [3]). The mass-metallicity rela- 
tions predicted by our best models with different numbers 
of bursts are shown in Fig. 1161 and they are consistent with 
the observed one considering the scatter of the data, espe- 



cially in the cases where the wind develops in the whole 
galactic mass range (i.e., n > 5). 

Considering that dwarf galaxies could be in different 
evolutionary stages and/or have different ages (measured 
from the beginning of star formation), we show that our 
best models - a series of uniform models with same Mi n f , 
wh^c, X mw , £) d but different n and suitable t - can well 
reproduce the spread in the observations. Therefore, our 
preferred galaxy formation scenario for these galaxies is 
the following: they should have accreted a lot of primor- 
dial gas at their early stages, and formed stars through sev- 
eral short star bursts (d ~ 0.3 Gyr for Afj n j = 10 9 Mq), 
with more massive galaxies suffering longer star formation 
bursts. However, it is likely that the gas escape from the 
potential well of the galaxy when enough energy from SN 
explosions is accumulated, and this wind should be metal- 
enhanced (u>h,Hc ~ 0.3/ The wind rate should be propor- 
tional to the gas mass at that time, and the wind efficiency 
should be X mw < 1. 

DLAs could be the progen itors of dwarf irre g ular g alax- 
ies, as already sugg ested by iMatteucci et al. I (|1997fl and 
ICalura et al.T(l20ol . In our best models, we adopted a 
fixed duration for each burst with the purpose of changing 
the parameters as less as possible. The evolutionary tracks 
of our best models pass through the DLA data, but could 
not explain the scatter. However, if we reduce the duration 
of the first SF burst to d = 0.01 ~ 0.1 Gyr, by taking model 
M9b3 as an example, the models can explain the scatter in 
the abundance ratios of DLAs much better (see Fig. [17] and 
Table [3] for details of these models) . Therefore, we confirm 
that the DLA systems could be progenitors of local dwarf 
galaxies. 

4.5. Continuous star formation 

By using a spectroph otometic model coupled with a chem- 
ical evolution model, iLegrandl (|2000l ) demonstrated that a 
continuous but very mild star formation rate (SFE as low 
as 10 -3 M Gyr _1 ) is able to reproduce the main proper- 
ties of IZw 18, one of the most metal-poor BCDs we know. 
Therefore, in this section we are showing model results ob- 
tained with continuous and mild SF. 

We run the same models as in ILegrand I (|200(J, one 
model with a mild continuous SF (e = 0.001 Gyr -1 ) only, 
and the other one with a mild continuous SF (e = 0.001 
Gyr -1 ) and a current burst (SFR = O.O23M yr" 1 , i.e., e = 
0.88 Gyr -1 if the observed M m = 2.6 x 1O 7 M , during the 
last 20 Myrs). The abundance ratios of different elements 
as predicted by these two models are shown in Fig. [T5J 
As expected, such a low continuous SF predicts a too low 
oxygen abundance which could not explain the majority 
of local BCDs. In addition, by examining the evolutionary 
tracks in the low metallicity range, one sees that DLAs 
cannot be the progenitors of such galaxies. 

Therefore, we developed other models with higher SFEs 
and different SF duration. Since both dlrrs and BCDs har- 
bor recent SF activities, we assumed that the SF is still 
going on at the present day, and therefore a short duration 
d implies a late starting time (to — d) of SF. 

We have examined the models with different SFHs by 
varying the SFE, duration, and infall timescale. The re- 
sults show that only models with different SFEs are able 
to explain the scatters in the abundance ratios. In Fig. \T§\ 
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Fig. 17. The evolutionary track of abundance ratios as predicted by our DLA models. All of these models assume same 
SFHs as model M9b3 except for a shorter duration of the first burst used. The red dotted, short-dash, long-dash and solid 
lines are for models DLAl to DLA4 whose durations of the first burst are d = 0.01,0.02,0.05,0.1 Gyr respectively. As 
a comparison, we also show the best model M9b2 of dwarf galaxy in black-dash-dot lines. All these models assume the 
same total infall mass Mj„/ = 10 9 Mq. Cyan open circles are BCDs and green open squares are DLAs, same as in Fig. [31 



models with different SFEs for the cases of no wind (left col- 
umn), normal wind (middle column), and metal-enhanced 
wind (wh.Ho = 0.3, right column) are shown. In these mod- 
els, to avoid reaching a too high metallicity at the present 
time, a duration shorter than the age of the universe is as- 
sumed for the highest SFE case. A short infall timescale is 
adopted (r = 1 Gyr) here, because it does not effect the 
evolutionary track too much in the continuous SF scenario, 
especially when the SF starts late. 

There are several conclusions that can be drawn from 
Fig. [191 

1. In the case of no wind, although models with different 
e could partly explain the scatters in the abundance 
ratios, they cannot reproduce the \x — Z relation. As we 
can see from the bottom left panel of Fig. [T9l all the 
evolutionary tracks are overlapping, the same as in the 
bursting SF scenario without galactic wind; 



2. When the wind is included in the model, the gas frac- 
tion decreases with time (the bottom middle and right 
panels of Fig. IT9"|) . We adopt relatively high wind effi- 
ciencies in the continuous SF scenario. When we com- 
pare the models with the same SFE in both normal 
wind and metal-enhanced wind cases, we can see that 
their gas fractions reach almost the same values at the 
present time since their gas loss rates are comparable 
(X w = Wn t jn e X mw = 0.6); however, their metals show 
different behaviors. The oxygen abundance decreases 
dramatically in the metal-enhanced wind case, thus ex- 
plaining the scatter in fi — Z relation much better. In the 
meantime, the metal enhanced wind model can fit the 
data relative to abundance ratios much better (upper 
and middle panels of right column in Fig. [Ti?)k 

3. As we have demonstrated before, galaxies with long con- 
tinuous but mild SF (e.g., e = 0.01, d = 13), which 
could be the case for some BCDs, as suggested by 
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Fig. 18. The evolutionary track of abundance ratios as predicted by the models of Leerand ( 2000h . A mild continuous 
SF with e = 0.001 Gyr -1 is shown in black-solid lines, whereas a mild continuous SF (e = 0.001 Gyr -1 ) combined with 
a current burst (SFR = O.O23M yr- 1 , i.e., e = 0.88 Gyr" 1 if M H i = 2.6 x 10 7 M o , during the last 20 Myrs) is shown in 
red-dash lines. The blue shade areas show the continuous SF wit h 1.45 x 10~ 3 < e < 3.85 x 10 -3 Gyr -1 which is derived 
from the observed SFR and A/hi data of IZw 18(see Table [Tl of iLeerand" 200Clf h Cyan open circles arc BCDs and green 
open squares are DLAs, same as in Fig. [3] 



the p revious work (e.g., iLegrandl 120001 : ILeerand et al. I 
l2000f ) . should not be the majority due to the too low 
predicted oxygen abundances. Their evolutionary tracks 
relative to abundance ratios also do not fit the data. In 
addition, DLAs cannot be the progenitors of these ob- 
jects when abundance ratios are examined, because of 
the very high predicted values of N/O, N/Fe and N/Si 
at low metallicity (sec Fig. IT5|) . at variance with the 
properties of DLAs. 
4. In the continuous SF scenario, the optimal models 
should have e = 0.05 — 0.2, a duration of SF shorter 
than 13 Gyr, and metal-enhanced winds should occur, a 
situation similar to a long starburst scenario. Moreover, 
DLAs can be the progenitors of these galaxies. 



5. Discussion and Conclusions 

We have discussed in detail the chemical evolution of late- 
type dwarf galaxies (dwarf irregular and blue compact 
galaxies) and used the most recent data as a comparison. 
We have taken into account the measured abundances of 
single elements (He,C, N, O, S, Si and Fe) as well as the gas 
masses. We have assumed that the late- type dwarf galaxies 
form by cold gas accretion and we run models for differ- 
ent accreted baryonic masses (10 8 ,10 9 and 1O 1O M0). We 
have tested both bursting and continuous star formation. 
We then have studied in detail the development of galactic 
winds by assuming a dark matter halo which is assumed 
10 times the amount of of the baryonic mass and feedback 
from SNe and stellar winds. Our main conclusions are: 

1. Galactic winds are necessary to reproduce the main 
properties of late-type dwarf galaxies and they should 
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Fig. 19. The evolutionary track as predicted by models with continuous star formation. From top to bottom panels are 
log(C/0) vs. 12+log(0/H), [O/Fe] vs. [Fe/H], and the mu — Z relation; from left to right columns are for the cases of 
no wind, normal wind, and metal-enhanced wind {wn,He = 0.3). Models with different SFEs are shown in each case: 
black-solid lines for e = 0.5 {d = 2), magenta- dash- dot lines for e = 0.2 (d = 5), blue-long- dash lines for e = 0.1 (d = 7), 
red-dash lines for e = 0.05 (d = 9), and black-dot lines for e = 0.02 (d = 13). The total infall masses are assumed to be 
Minf = 10 9 A/ Q in all these models. The observational data are the same as in Fig. [5] and Fig. [31 BCDs, dlrrs, and DLAs 
arc plotted in cyan open circles, cyan open triangles and green open squares respectively. 



be metal-enhanced, namely metals should be carried 
away preferen tially, in agreement with previous dynam- 
ical w ork (e.g. iMac Low fc Ferrara | [l999: Re cchi et al. I 
l200lh . The rate of gas loss was assumed to be propor- 
tional to the amount of gas present at the time of the 
wind, which is equivalent to say that is proportional to 
the SFR, in agreem ent with previous papers and ob- 
servational eviden ce (jMartin Il20"05l: iRupke et aT~ll2005t 
IChen et al. Il2010l ). The wind efficiency X mw should be 
relatively low (A mul ~ 0.8 in bursting SF scenario, and 
X mw ~ 2 in continuous SF scenario) compared to what 
is assumed for dwarf spheroidals (~ 6 to 15) where 
the wind should carry away a ll the residual gas (e.g. 
lLanfranchi fc Matteucci II2003I ). 
2. Both bursting and continuous SF scenarios have been 
examined. In the case of bursting SF, the number of 



bursts could be no more than ~ 10 and the star forma- 
tion efficiency could be ~ 0.5 Gyr -1 . Galaxies with a 
long continuous but mild SF (e < 0.02 Gyr -1 , d ~ 13 
Gyr) should not be the majority, whereas galaxies with 
higher SFE (0.05 Gyr" 1 < e <0.2 Gyr" 1 ) and shorter 
SF duration (5 Gyr < d < 9 Gyr) are still acceptable. 

3. Models with a different number of bursts and/or dif- 
ferent star formation efficiency and/or different burst 
duration, coupled with metal enhanced winds, can re- 
produce at best the spread observed in the abundances 
versus fractionary mass of gas and abundance ratios ver- 
sus abundances. Normal winds, where all the gas and 
metals are lost at the same rate, should be rejected since 
they subtract too much gas. 

4. We studied the M — Z relation for late- type dwarf galax- 
ies and we showed that in order to reproduce such a re- 
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lation one should assume that different galaxies suffered 
an increasing amount of star formation (efficiency, num- 
ber of bursts, duration of bursts) with galactic mass, 
even without galactic winds. However, metal-enhanced 
wind are necessary to reproduce all the other features. 
In the case of the M — Z relation, a metal-enhanced 
wind efficiency increasing with galactic mass can very 
well reproduce the data, leaving all the other param- 
eters to be the same irrespective of the galactic mass. 
On the other hand, normal winds occurring at a dif- 
ferent efficiency in different galaxies cannot reproduce 
the M — Z relation unless the primordial gas is infalling 
continuously (i.e., a long infall timescale), or the SFE 
or the number of bursts increase with galactic mass. 

5. A comparison of our best model predictions with data 
for DLAs has shown that these objects can well be the 
progenitors of local dlrrs and BCDs, in agreement with 
previous papers. 

6. In this work, our models can reproduce the chemical 
properties of both dlrrs and BCDs. To distinguish be- 
tween these two types of galaxies, the photometric and 
spectral information should also be taken into account. 
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